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a Recently, an accurate numerical study of the crossover from asymptotic (Ising-like) critical behavior to classical 
(mean-field) behavior has been performed both for two-dimensional Jl],|| and three-dimensional |j| Ising systems in 
+Jj \ zero field on either side of the critical temperature with a variety of interaction ranges. It is the objective of the present 
work to analyze these numerical results within the framework of a crossover theory that is based on renormalization- 
group matching and that has already successfully been applied to the description of crossover in several experimental 
systems [|,|. 

Qualitatively the crossover is ruled by the parameter t/G where t — (T—T c )/T c is the reduced temperature distance 
to the critical temperature T c and G the Ginzburg number [[|. The Ginzburg number depends on the normalized 
interaction range R as 
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G = G R- 2d/(4 - d) , (1) 



where d is the dimensionality of space and Go a constant. Hence, for d = 3 the crossover occurs as a function of tR 6 . 
Asymptotic critical behavior takes place for tR 6 <C 1 and classical behavior is expected for tR 6 3> 1. In real fluids the 
crossover is never completed in the critical domain (where t <C 1), since the range of interaction is of the same order 
, of magnitude as the distance between molecules (R ~ 1) |Q. A new Monte-Carlo algorithm, developed by Luijten 
ly-j ' and Blote 0, offers the advantage that the ratio t/G can be tuned over more than eight orders of magnitude allowing 
, one to cover the full crossover region in three-dimensional spin models || . 

A sensitive description of crossover behavior is obtained from an analysis of the effective critical exponent of the 
\ susceptibility (the third derivative of the free energy), defined as 

7 ± =-dlnx/dln\t\, (2) 

S where the scaled susceptibility x — kBT c (R)(dm/dh)T, fee the Boltzmann constant, m the order parameter, h the 
ordering field, and where the "+" sign applies for T > T c , and the "— " sign for T < T c . As is seen from Figs. 1 and 2, 
, the variation of 7^ reproduces the Ising asymptotic critical behavior (7^ ~ 1.24) at tR 6 <C 1 as well as the mean-field 
asymptote (7^ = 7mf = 1) at tR 6 3> 1. Apparently, all data would seem to collapse onto a universal function of the 
reduced variable tR 6 as predicted by a field-theoretical treatment |||| and by the e-expansion pQ] . However, as was 
noted in Ref. || , a more careful look at the data reveals a remarkable discrepancy between the theoretical calculations 
]l2| and the simulation results. Namely, the shape of the crossover is sharper than predicted by the theory |jil,[l2 |, 



especially for short ranges of interaction. We will show that this discrepancy is related to the findings of Refs. 



where it was shown that there is a fundamental problem in describing the crossover of 7^ by a universal function 
^ ■ which contains only a single crossover parameter G oc R~ e . 

In zero-ordering field above T c the susceptibility asymptotically close to the critical point behaves as 

x = r <- 7 (i + rdi A * + r 2 < 2A * + ai t + ...), (3) 
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where 7 = 1.239 ± 0.002 (see, e.g., Refs. HQ and references therein) and A s = 0.504 ± 0.008 |y| 

are universal 

Ising critical exponents, and where Tq, T%, T2, and a\ are system-dependent amplitudes. Expansion^ is called the 
Wegner series [ JL6| . 

In a universal single-parameter crossover theory [^|-^o|, the Ginzburg number is responsible both for the range of 
validity of the mean-field approximation and for the convergence of the Wegner series (^). However, it is known 
fl7| |l9| , that the sign of the first Wegner correction amplitude Ti depends on the difference u — u* , where u is the 
scaled coupling constant and u* — 0.472 is the universal coupling constant at the Ising fixed point ^0|. Moreover, 
Liu and Fisher fl8f| concluded that the three-dimensional nearest-neighbor Ising model has a negative leading Wegner 
correction amplitude I\, so that 7^ asymptotically approaches 7 ~ 1.24 from above. Therefore, since the coupling 
constant itself depends on the interaction range, the shape of 7^ cannot be represented by a universal function of 
the Ginzburg number, since G is not proportional to the difference u — u* . 

In this paper we therefore present an analysis of the numerical data for 7^ |3| in terms of a crossover model 
based on renormalization-group matching for the free-energy density |l7|Jl9| , |2l| . This model contains two crossover 
parameters u = u/u* and A (a dimensionless cutoff wave number), and two rescaled amplitudes Ct and c p related to 
the coefficients of the local density of the classical Landau-Ginzburg free energy AA: 
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with t = (T - T c )/T, M = c p (p = (a /c t ) 1/2 (p, a Q 



Uq = u*uAcp, 
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and V = fg^ 3 V. The average 



molecular volume vq and the prefactor vo/ksT are introduced to make the free-energy density and all the coefficients 
dimensionless. The inverse crossover susceptibility x^ 1 = (d 2 AA/dM 2 ^J , where AA is the crossover (renormalized) 
free-energy density, in zero field above T c reads H 



X- 1 =c 2 p c t TY^/^(l + y) 



(5) 



with 



2A S 



v (l-u)Y 



1 - (1 - u)Y 



2v- 1 



(6) 



where v ~ 0.630 15 2^] is the critical exponent of the asymptotic power law for the correlation length £ 0]. Note that 
X -1 = (r c /T)x _ ^ind the relation between j cS = ~d\nx/d\n \t\ and 7^, given by Eq. (||), is 7^ = j eS + (1 - %s)r, 
both above and below the critical temperature. The crossover function Y is defined by 



1 - (1 - u)Y 
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and is to be found numerically. The parameter k in Eq. (^) is inversely proportional to the fluctuation-induced 
portion of the correlation length and serves as a measure of the distance to the critical point. In zero field above T c 
the expression for n 2 reads: 



(8) 



We modified the original expression for n 2 , given by Eq. (3) in [Q, by introducing the non-asymptotic factor T/T c in 
Eq. (8) so that k 2 becomes infinite at T — * 00 [^3|. Asymptotically close to the critical point (A/k^$> 1), the following 
expression is obtained for the first correction amplitude Ti in Eq. (0): 
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where g\ ~ 0.62 is a universal constant fl2l| . 

In the approximation of an infinite cutoff A — > 00, which physically means neglecting the discrete structure of 
matter, u = uoc 2 / (u* Aoq) — > and the two crossover parameters u and A in the crossover equations collapse into a 
single one, uA, which is related to the Ginzburg number G by Q 
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G = g^-^=gv , " " , (10) 



where 50 — 0.028 is a universal constant [^l| and £0 = v o^ c 7 1 ^ 2 = (co/ao) 1 ^ 2 is the mean-field amplitude of the power 
law for the correlation length. Note that the Ginzburg number does not depend explicitly on the cutoff A or on u. 
This single-parameter crossover, i.e., the crossover for u = 0, is universal and is indicated in Fig. 1 by a dashed-dotted 
curve. This simplified description of the crossover is equivalent to the results of Bagnuls and Bervillier || and of 
Belyakov and Kiselev E0 |. 

In the simulations |J]7 each spin interacts equally with its z neighbors lying within a distance R m on a three- 
dimensional cubic lattice. The effective range of interaction R is then defined as R 2 — z X^-=« \ r i ~ r j\ with 
l r i — r j| < Rm [0- We have approximated the relation between R and R m by R 2 = §.^(1 + as indicated 

in the insert in Fig. 3. In order to compare the numerical results to the theoretical prediction Eq. (|5|), we need the 
range dependence of the parameters c t and u. Indeed, the asymptotic R dependence of u follows directly from simple 
scaling arguments [pi , u = uoR~ 4 , and Ct varies as its square root, c t = c t oR~ 2 . For a three-dimensional simple cubic 
lattice, A = ir [[L8||24|] , and we obtain for the Ginzburg number 

G = G iT 6 = 0.28(^/4,)^ = 0.28(u 2 /c t0 )R- 6 . (11) 

The nonuniversal parameters Cto and uq have to be determined from a least-squares fit to the numerical data for 7^, 
which yielded c t o — 1.72 and uq = 1.22 and hence Go ~ 0.24. The solid lines in Fig. 1 indicate the corresponding 
theoretical curves. It should be noted that these curves are calculated for each value of R m separately; the piecewise 
continuous character of this description directly reflects the fact that the crossover cannot be described by a universal 
single-parameter function. Indeed, Fig. 1 also shows two attempts to describe the data in terms of such a function. 
The dash-dotted line corresponds to the limit u — > 0, whereas the dotted curve corresponds to u$ = 1.22 and A = ir 
(a continuation of the theoretical curve for R = 1). We see that the actual crossover lies between these two bounding 
curves, with u ~ for large R and u ~ 1.2 for R = 1. Thus, it is clearly seen that without including the R dependence 
of u it is impossible to describe data for short interaction ranges R 2 U < 5. The dependence of u on R is shown in 
Fig. 3. The two adjustable parameters Cto and uq are strongly correlated and if one of them is fixed at a predicted 
value, the quality of the description remains the same. We hence tried to fit the data while keeping cto fixed at the 
theoretically predicted value Cto = 2d = 6 |2^,^]. In this case a fit of the same quality is obtained with uq = 1.22, 
provided that A ~ 2ir. The value of Go ~ 0.24 then remains unchanged. 

To describe the data below the critical temperature, a connection between M and t in zero field is to be found from 
the condition (dAA/dM) r = 0. The relation between M and r appears to be implicit and x as a function of r cannot 
be expressed in an explicit form either. Of course, the parameters Cto and uq should be the same as for T > T c and we 
hence kept them fixed at the above-mentioned values. However, the parameter Go appearing in Eq. ( pi] ) will take a 
different value. We took this into account by introducing a factor Gq /Gq into the temperature scale: t — > t ■ (Gq / Gq ) ■ 
Figure 2 shows the results for T < T c , where the factor Gq /Gq was included as an adjustable parameter. Our estimate 
Gq /Gq — 2.58 must be compared with the theoretical result Gq /Gq = 3.125 p7|| . Interestingly, 7^ clearly shows a 
minimum around |i|i? 6 ~ 10 2 . This corroborates the nonmonotonic character of the crossover of 7^, earlier observed 
for the two-dimensional Ising lattice |2j, where the effect is much more pronounced. We note that already in Ref. |2S[] 
a field-theoretic calculation of the crossover in the low-temperature re gim e has been given (in the limit u — * 0), but 
only recently this has been extended to cover the full crossover region p9[ . Actually, also here a nonmonotonicity in 
7^ has been observed. 

In summary, we remark that although in general the theory contains two crossover parameters u and A, only one 
parameter (u) changes with the range of interaction. However, this does not mean that the crossover is a universal 
function of tR 6 . Indeed, the effective range of interaction R affects the behavior of 7^ in a twofold way: through the 
Ginzburg number, which is proportional to c)?, and through the first Wegner correction, with an amplitude Ti that 
is proportional to (1 — u) [Eq. (0)]. Hence, there is no way to describe the data for short interaction ranges without 
allowing for u to become larger than unity and correspondingly Ti to change its sign between i? m = 2 and R m = 1 
as indicated in Fig. 3. In previous publications we have shown that Eq. (||), derived from renormalization-group 
matching, gives an excellent representation of the experimentally observed crossover behavior in simple and complex 



fluids Mp 30 1 . From the evidence presented in this paper, we conclude that the same crossover model also yields a 



quantitative description of the crossover critical behavior of a three-dimensional Ising lattice. 
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FIG. 1. The effective susceptibility exponent 7^ above T c . The symbols indicate numerical simulation data ||. The solid 
curves represent values calculated from Eq. Q). The dashed-dotted curve corresponds to the limit u — > 0. The dotted curve 
is a continuation of the crossover curve for u = 1.22. For clarity, the error bars have been omitted; they are all of the order 
of 0.004. 
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FIG. 2. The effective susceptibility exponent 7^ below T c . The symbols indicate numerical simulation data 
curves represent values calculated from the renormalization-group matching crossover model. 
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FIG. 3. Dependence of the normalized coupling constant u on the normalized interaction range R. Note that u becomes 
larger than unity for very short interaction ranges. Insert: Effective range of interaction R (open circles) plotted as a function 
of Rm- The solid line corresponds to the approximation mentioned in the text and the dashed line represents the asymptotic 
behavior for large R. 
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